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Abstract 

We  evaluated  the  performance  of  system  combining  a  solid  oxide  fuel  cell  (SOFC)  stack  and  a  polymer  electrolyte  fuel  cell  (PEFC)  stack 
by  a  numerical  simulation.  We  assume  that  tubular-type  SOFCs  are  used  in  the  SOFC  stack.  The  electrical  efficiency  of  the  SOFC-PEFC 
system  increases  with  increasing  oxygen  utilization  rate  in  the  SOFC  stack.  This  is  because  the  amount  of  exhaust  heat  of  the  SOFC  stack 
used  to  raise  the  temperature  of  air  supplied  to  it  decreases  as  its  oxygen  utilization  rate  increases  and  because  that  used  effectively  as 
the  reaction  heat  of  the  steam  reforming  reaction  of  methane  in  the  stack  reformer  increases.  The  electrical  efficiency  of  the  SOFC-PEFC 
system  at  190  kW  ac  is  59%  (LHV),  which  is  equal  to  that  of  the  SOFC-gas  turbine  combined  system  at  1014  kW  ac. 
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1.  Introduction 

The  electrical  efficiency  of  a  power  generation  system  us¬ 
ing  a  solid  oxide  fuel  cell  (SOFC)  stack  is  expected  to  be 
very  high  [1].  An  SOFC  cogeneration  system  that  generates 
electricity  and  heat  used  for  heating  water  has  achieved  46% 
(LHV)  electrical  efficiency  at  109  kW  ac  [2].  However,  it 
has  been  pointed  out  that  an  SOFC  system  that  generates 
electricity  in  the  SOFC  stack  and  other  generating  equip¬ 
ment  can  provide  higher  electrical  efficiency  than  the  SOFC 
cogeneration  system  [1]. 

SOFC  systems  that  generate  electricity  in  an  SOFC 
stack  and  gas  turbine  (GT)  are  receiving  attention  [3,4]. 
At  19.5  MW  ac,  the  electrical  efficiency  of  a  SOFC-GT 
system  is  expected  to  reach  70%  [3].  However,  SOFC-GT 
systems  are  lauder  than  SOFC  cogeneration  systems,  and 
the  cost  of  the  auxiliary  equipment  is  higher.  In  addition, 
the  energy  conversion  efficiency  of  the  GT  decreases  as  the 
output  decreases,  so  at  220  kW  ac  the  electrical  efficiency 
of  a  SOFC-GT  system  using  a  50kW-class  GT  is  expected 
to  be  57%  [4],  which  is  far  lower  than  that  at  19.5  MW  ac. 
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Systems  that  generate  electricity  in  the  SOFC  stack  and  a 
polymer  electrolyte  fuel  cell  (PEFC)  stack  are  also  attracting 
attention  [5,6].  The  SOFC-PEFC  system  can  provide  higher 
efficiency  than  SOFC  cogeneration  systems  [5,6]  like  the 
SOFC-GT  system.  The  SOFC-PEFC  system  is  quieter  than 
the  SOFC-GT  system,  and  the  cost  of  the  auxiliary  equip¬ 
ment  is  lower.  The  SOFC-PEFC  systems  can  provide  higher 
efficiency  than  the  SOFC-GT  systems  when  the  output  is 
comparatively  small,  because  the  efficiency  of  the  PEFC 
stack  remains  almost  constant  even  as  the  output  decreases. 
However,  the  performance  of  the  SOFC-PEFC  system  has 
not  been  clarified  yet. 

We  quantitatively  simulated  and  clarified  the  performance 
of  a  200kW-class  SOFC-PEFC  system  with  a  lOOkW-class 
tubular  SOFC  stack. 


2.  Simulation  models 

2.7.  SOFC  stack 

The  SOFC  stack  consists  of  1152  stack  units,  each  of 
which  consists  of  a  tubular  SOFC,  stack  reformer,  pre¬ 
reformer,  and  combustor  as  shown  in  Fig.  1.  The  tubular 
SOFC  is  composed  of  the  anode,  cathode,  and  solid  oxide 
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Nomenclature 

A 

heat  transfer  area  (m2) 

C 

isopiestic  specific  heat  (J/mol  K) 

d 

thickness  (m) 

E 

electromotive  force  (V) 

F 

Faraday’s  constant  (C/mol) 

h 

heat  transfer  coefficient  (W/m2  K) 

J 

current  density  (A/m2) 

K 

equilibrium  constant 

l 

length  (m) 

L0 

energy  loss  ratio  to  gross  ac  output 

M 

molar  flow  rate  (mol/s) 

Ma 

molar  flow  rate  of  anode  exhaust  gas 
from  tubular  SOFC  at  minimum  oxygen 
utilization  rate  to  keep  Tm  at  setting 
temperature  (mol/s) 

N 

number  of  tubular  SOFCs  in  SOFC 

stack 

Nu 

non-dimensional  heat  transfer  coefficient 

P 

partial  pressure  (Mpa) 

Q 

radiation  loss  from  100  kW  class  SOFC 
stack  made  by  SWPC  (W) 

R 

gas  constant  (J/mol  K) 

Rb 

base  resistance  (Q  m2) 

Rc 

contact  resistance  between  each  tubular 
SOFC  (Qm2) 

Re 

summation  of  anode,  cathode,  and 
electrolyte  resistance  of  tubular 

SOFC  (Qm2) 

s 

cross-sectional  area  (m2) 

T 

temperature  (K) 

Th 

temperature  for  definition  of  enthalpy 
change  (K) 

Tm 

maximum  temperature  in  tubular 

SOFC  (K) 

Uf 

fuel  utilization  rate  (%) 

U0 

oxygen  utilization  rate  (%) 

V 

cell  voltage  (V) 

Vv 

volume  of  pre-reformer  (m3) 

w 

ac  output  (W) 

X 

coordinate  along  gas  flow  (m) 

Greek  letters 

AH 

enthalpy  change  at  Tbase  (J/mol) 

xj/ 

dc/ac  conversion  efficiency 

8 

internal  diameter  (m) 

£ 

activation  energy  (J/mol) 

0 

ohmic  overpotential  of  tubular  SOFC  (V) 

0 

electrical  efficiency  (%) 

<P 

external  diameter  (m) 

K 

frequency  factor 
(mol/(m3/cat)s  MPa09845) 

X 

heat  conductivity  (W/m  K) 

Subscript 

A 

air 

AEG 

anode  exhaust  gas 

ac 

alternating  current 

ALM 

alumina  tube 

AMB 

ambient 

ANO 

anode  of  tubular  SOFC 

APE 

anode  exhaust  gas  for  PEFC  stack 

AUX 

auxiliary  machine 

AVE 

average 

BLOW 

air  blower 

CATH 

cathode  of  tubular  SOFC 

CEG 

combustion  exhaust  gas 

CELL 

anode,  cathode,  and  electrolyte  assembled 

solid 

COM 

combustor 

dc 

direct  current 

ELE 

electrolyte 

F 

fuel 

FCHA 

fuel  channel 

HE 

heat  exchanger 

IN 

inlet 

INS 

inside 

OTHER 

other  than  air  blower 

OUT 

outlet 

OUTS 

outside 

OXI 

oxidation 

PE 

PEFC 

PEstack 

PEFC  stack 

PRE 

pre-reformer 

REC 

recycled  gas  for  steam  reforming 

SHIF 

shift  reaction 

SO 

tubular  SOFC 

SOstack 

SOFC  stack 

SR 

steam  reforming  reaction 

SRE 

stack  reformer 

SYS 

SOFC-PEFC  system 

Wall 

wall 

electrolyte.  We  refer  the  anode,  cathode,  and  electrolyte  as¬ 
sembled  solid  as  the  cell  in  this  paper.  The  tubular  SOFCs 
are  electrically  connected  to  each  other  through  Ni  felt  in 
the  SOFC  stack  [7].  The  methane  fed  to  the  stack  unit  is 
mixed  with  the  recycle  gas  for  steam  reforming  and  sup¬ 
plied  to  the  pre-reformer.  Some  of  the  methane  is  converted 
to  hydrogen  and  carbon  monoxide  by  the  steam  reform¬ 
ing  reaction  in  the  pre-reformer.  The  exhaust  gas  from 
the  pre-reformer  is  supplied  to  the  stack  reformer,  where 
the  residual  methane  is  converted  to  hydrogen  and  carbon 
monoxide  by  the  steam  reforming  reaction  using  the  ex¬ 
haust  heat  of  the  tubular  SOFC  as  the  reaction  heat.  The 
exhaust  gas  from  the  stack  reformer  is  supplied  to  the  an- 
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AC  output 

• — 


,  auxiliary  machine  for  SOFC  stack 


gross  AC  outputl  inverter 
DC  output 


methane 


anode  exhaust  gas  for  PEFC  stack 
from  SOFC  stack 


SOFC® 

/ 

shift 

CO  selective 

stack 

4 

converter 

► 

oxidizer 

heated 

air 


exhaust 

combustion  gas 
from  SOFC  stack 

heat 

exchanger 


condenser 


exhaust  gas 


water 


<s 

> 

V 

cathode^ 

PEFC 
stack  t 

exhaust  gas 


anode  exhaust  gas  from  PEFC  stack 


DC  output 
inverter 

gross  AC  output 
AC  output  • 


E 

0 


auxiliary  machine  for  PEFC  stack 


ode  of  the  tubular  SOFC  and  used  for  power  generation.  A 
part  of  the  anode  exhaust  gas  is  discharged  from  the  stack 
unit  as  the  anode  exhaust  gas  for  the  PEFC  stack.  All  anode 
exhaust  gases  for  the  PEFC  stack  from  the  stack  units  are 
gathered  and  supplied  to  the  shift  converter  as  the  anode 
exhaust  gas  for  the  PEFC  stack  from  the  SOFC  stack  (see 
Fig.  2).  Another  part  of  the  anode  exhaust  gas  from  the 
tubular  SOFC  is  recycled  to  the  pre-reformer  as  the  recycle 
gas  for  steam  reforming.  We  chose  a  molar  flow  rate  of 
the  recycle  gas  for  steam  reforming  so  that  the  steam  to 
methane  molar  ratio  (S/C  ratio)  in  the  gas  supplied  to  the 
pre-reformer  is  3.0.  The  rest  of  the  anode  exhaust  gas  from 
the  tubular  SOFC  is  fed  to  the  combustor  together  with  the 
anode  exhaust  gas  from  the  PEFC  stack  and  the  cathode 
exhaust  gas  from  the  tubular  SOFC.  The  methane,  hydro¬ 
gen,  and  carbon  monoxide  in  the  two  anode  exhaust  gases 
completely  react  with  the  oxygen  in  the  cathode  exhaust  gas 
from  the  tubular  SOFC  in  the  combustor.  All  combustion 
exhaust  gases  from  the  combustors  of  the  stack  units  are 
gathered  and  fed  to  the  heat  exchanger  as  the  combustion 
exhaust  gas  from  the  SOFC  stack  (see  Fig.  2).  The  air  is 
fed  to  the  alumina  tube  installed  inside  the  tubular  SOFC, 
and  the  temperature  of  the  air  is  raised  by  using  the  ex¬ 
haust  heat  of  the  tubular  SOFC.  The  heated  air  is  supplied 
to  the  cathode  of  the  tubular  SOFC  and  used  for  power 
generation. 

We  used  the  following  assumptions  for  the  simulation  of 
the  SOFC  stack. 

(1)  The  methane,  air,  and  anode  exhaust  gas  from  the  PEFC 
stack  fed  to  the  SOFC  stack  are  distributed  to  every 
stack  unit  equally. 


Fig.  2.  SOFC-PEFC  system. 

(2)  Both  the  steam  reforming  reaction  of  methane, 

CH4  +  H2O^CO  +  3H2,  (1) 

and  the  shift  reaction, 

C0  +  H20*>C02  +  H2,  (2) 

occur  in  the  pre-reformer  and  stack  reformer.  The  shift 
reaction  is  in  thermodynamic  equilibrium  in  both  the 
pre-reformer  and  stack  reformer  [8]. 

(3)  The  oxidation  of  hydrogen, 

H2  +  \02  H20,  (3) 

the  oxidation  of  carbon  monoxide, 

CO  +  ±02  **  C02,  (4) 

and  the  oxidation  of  methane, 

CH4  +  202  2H20  +  C02,  (5) 

occur  in  the  tubular  SOFC  as  cell  reactions.  These  oxi¬ 
dations  are  in  thermodynamic  equilibrium  at  the  anode 
side  of  the  tubular  SOFC  [9]. 

(4)  The  temperature  profiles  of  the  anode,  cathode,  and  elec¬ 
trolyte  are  the  same. 

(5)  The  voltage  drop  in  a  tubular  SOFC  is  caused  by  over¬ 
potential  and  contact  resistance  between  tubular  SOFCs. 
The  overpotential  consists  of  only  ohmic  overpoten¬ 
tial  because  the  ohmic  overpotential  is  the  largest  one 
among  the  ohmic,  activation,  and  concentration  overpo- 
tentails  in  a  tubular  SOFC  made  by  Siemens  Westing- 
house  Power  Corporation  (SWPC)  [10].  The  logarithm 
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of  the  total  resistance  of  the  tubular  SOFC,  which  is  the 
summation  of  the  anode,  cathode,  and  electrolyte  resis¬ 
tances,  is  expressed  by  a  linear  function  of  the  inverse  of 
the  cell  temperature  because  the  logarithm  of  the  cath¬ 
ode  resistance,  which  is  the  largest  one  among  three  in 
the  tubular  SOFC  made  by  SWPC  [1 1],  is  expressed  by 
a  linear  function  of  the  inverse  of  the  cathode  tempera¬ 
ture  [12]. 

(6)  Heat  is  only  radiated  from  the  combustor. 

(7)  All  tubular  SOFCs  have  the  same  characteristics  in  the 
SOFC  stack. 

(8)  The  average  current  density  and  oxygen  utilization  rate 
in  the  SOFC  stack  are  the  same  as  those  in  the  tubular 
SOFC. 

2.2.  SOFC-PEFC  system 

The  configuration  of  the  SOFC-PEFC  system  is  shown 
in  Fig.  2.  As  mentioned  above,  the  anode  exhaust  gas  for 
the  PEFC  stack  from  the  SOFC  stack  is  supplied  to  the  shift 
converter,  and  the  exhaust  gas  from  the  shift  converter  is 
supplied  to  the  CO  selective  oxidizer.  The  carbon  monox¬ 
ide  in  the  anode  exhaust  gas  for  the  PEFC  stack  from  the 
SOFC  stack  is  completely  converted  to  carbon  dioxide  when 
the  gas  passes  through  the  shift  converter  and  the  CO  se¬ 
lective  oxidizer.  The  exhaust  gas  from  the  CO  selective  ox¬ 
idizer  is  supplied  to  the  condenser,  where  the  steam  is  cap¬ 
tured  as  water  by  decreasing  the  temperature  of  the  ex¬ 
haust  gas  from  the  CO  selective  oxidizer.  The  exhaust  gas 
from  the  condenser  is  supplied  to  the  PEFC  stack  and  hy¬ 
drogen  in  the  exhaust  gas  from  the  condenser  is  used  for 
power  generation.  As  mentioned  above,  the  anode  exhaust 
gas  from  the  PEFC  stack  is  distributed  to  each  stack  unit 
equally  and  burnt  in  the  combustor  of  the  stack  unit.  The 
combustion  exhaust  gas  from  the  SOFC  stack  supplied  to 
the  heat  exchanger  is  used  for  raising  the  temperature  of  the 
air  supplied  to  the  SOFC  stack.  The  dc  output  of  the  SOFC 
stack  and  that  of  the  PEFC  stack  are  inverter-converted 
ac  outputs.  We  call  them  gross  ac  outputs  in  this  paper. 
Some  of  the  gross  ac  outputs  are  used  for  the  auxiliary 
machine  for  the  SOFC  stack  and  the  one  for  the  PEFC 
stack. 

We  used  the  following  assumption  for  the  simulation  of 
the  SOFC-PEFC  system. 

(1)  Oxidation  of  hydrogen  expressed  as  Eq.  (3)  is  the  only 
cell  reaction  in  the  PEFC  stack  [13]. 

(2)  Ninety  nine  percent  of  the  carbon  monoxide  in  the  an¬ 
ode  exhaust  gas  from  the  SOFC  stack  is  converted  to 
carbon  dioxide  according  to  the  shift  reaction  expressed 
as  Eq.  (2)  in  the  shift  converter. 

(3)  The  carbon  monoxide  in  the  gas  supplied  to  the  CO  se¬ 

lective  oxidizer  is  completely  oxidized  to  carbon  diox¬ 

ide.  The  only  reaction  that  occurs  in  the  CO  selective 

oxidizer  is  the  oxidation  of  carbon  monoxide. 


3.  Fundamental  equations 

3.1.  SOFC  stack 

The  molar  flow  rate  of  methane  fed  to  the  pre-reformer, 
^ch4  -pre-in  ,  the  molar  flow  rate  of  air  supplied  to  the  alu¬ 
mina  tube,  Ma-alm-in,  the  dc  output  of  the  tubular  SOFC, 
Wso,  the  molar  flow  rate  of  the  anode  exhaust  gas  for  the 
PEFC  stack,  Mape,  and  the  molar  flow  rate  of  the  combus¬ 
tion  exhaust  gas,  Mceg,  are  calculated  using  the  average 
current  density  of  the  tubular  SOFC,  /so-ave,  the  oxygen 
utilization  rate  of  the  SOFC  stack,  £/0S0stack,  which  is  the 
same  as  that  of  the  tubular  SOFC,  and  the  maximum  tem¬ 
perature  of  the  tubular  SOFC,  Tm,  based  on  Eqs.  (6)-(47). 
The  /so-ave  and  Tm  are  2000  A/m1 2 3 * * *  [2]  and  1273  K  [14], 
respectively.  The  constants  used  in  the  simulation  are  listed 
in  Table  1. 

The  temperature  of  the  mixed  gas  supplied  to  the 
pre-reformer,  7f-pre-in,  is  calculated  by 

C  ch4  M ch4 -pre-in  (  7f-pre-in  —  Tamb) 

=  (Tf-SO-OUT  -  7e-pre-In)£  CjAfj-REC,  (6) 

i 

where  i  represent  methane,  hydrogen,  steam,  carbon  monox¬ 
ide,  and  carbon  dioxide. 

The  molar  flow  rate  of  methane  at  the  pre-reformer  outlet, 
4^ch4 -pre-out,  is  calculated  by  [15] 

4^ch4-pre-in  —  4^ch4-pre-out 

=  kSr  Vp  exp  ( —  —  )  Pch84-pre-in  •  (7> 

\F1  F-PRE-IN/ 

The  relationship  between  the  equilibrium  constant  of  shift 
reaction,  Ashif,  and  the  partial  pressure  of  each  gas  compo¬ 
nent  in  the  pre-reformer  is  expressed  as 

v  frr  v  PC02 -PRE-OUT  /?H2 -PRE-OUT 

F  SHIF  ( 1 F-PRE-OUT )  =  - •  (o) 

PCO-PRE-OUT  PH20-PRE-0UT 

The  mass  and  energy  balance  in  the  pre-reformer  is  ex¬ 
pressed  as 

(Mh2o-pre-in  —  4^h2o-pre-out) 

-  (4^ch4-pre-in  -  Mcra -pre-out) 

+  (Mco2-PRE-IN  4^co2-pre-out)  =  0,  (9) 

(Mh2-PRE-IN  -  4^h2-PRE-OUt) 

+  3  (Mch4 -PRE-IN  —  4^cH4-PRE-OUt) 

(Mco2-PRE-IN  4^co2-pre-out)  =  0,  (10) 

(Mco-PRE-IN  -  4^co-pre-out) 

+  (4^CH4-PRE-IN  —  3Tch4 -PRE-OUT ) 

+  (Mco2-PRE-IN  -  3Tco2 -PRE-OUT )  =  0, 


(11) 
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Table  1 

List  of  constants 


Constant 

Value 

ACom  (m2) 

2.44  x  10-3 

Ahe  (m2) 

10.7 

CCh4  (J/molK) 

61.9 

Cc0  (J/molK) 

31.8 

CCo2  (J/molK) 

50 

Ch2  (J/molK) 

29.9 

Ch2o  (J/molK) 

44 

Cn2  (J/molK) 

31.4 

Co2  (J/molK) 

33.3 

^ALM  (m) 

1  x  10“3 

<^ANO  (m) 

00 

1 

o 

X 

^CATH  (m) 

2.2  x  10“3  [18] 

Jele  (m) 

4  x  10-5  [18] 

^Wall  (m) 

1  x  10“3 

F  (C/mol) 

96484 

hcou  (w/m2  K) 

40  [19] 

hH E  (W/m2  K) 

40  [19] 

/cell  (m) 

1.2  [2] 

/fcha  (m) 

5.2  x  10“3  [2] 

ToAUX-Pestack 

0.93  [20] 

f-oBLOW 

0.04  [3] 

Toother 

0.03  [3] 

N 

1152  [2] 

N uA-ALM-INS 

4.36  [16] 

^VuA-ALM-OUTS 

9.98  [16] 

A^uA-CATH 

7.19  [16] 

A^uF 

8.23  [16] 

Q  (W) 

1.8  x  104  [2] 

Rr  (^m2) 

1 

Rc  (f2m2) 

3.5  x  10-5  [2,18] 

As  re  (m2) 

1.15  x  10-4  [2] 

Tamb  (K) 

298 

Th  (K) 

1173 

Vp  (m3) 

1.08  x  10-5 

<$alm  (m) 

6.8  x  10-3 

^so  (m) 

1.76  x  10“2 

A//ch4-oxi  (J/mol) 

-8.03  x  105 

A/Zco-oxi  (J/mol) 

-2.82  x  105 

A//h2-oxi  (J/mol) 

-2.49  x  105 

AZ/shif  (J/mol) 

-3.22  x  104 

AH$r  (J/mol) 

2.26  x  105 

ssr  (J/mol) 

2.28  x  104  [15] 

<PALM  (m) 

8.8  x  10“3 

<Pso  (m) 

2.2  x  10“2  [18] 

ksr  (mol/(m3/cat)sMPa0-9845) 

1.65  x  104  [15] 

Aa  (W/mK) 

0.078 

7-alm  (W/mK) 

2.6  [21] 

A-ano  (W/mK) 

6  [21] 

Tcath  (W/mK) 

11  [21] 

A-ele  (W/mK) 

2.7  [21] 

Af  (W/mK) 

0.51 

TSre  (W/mK) 

6.28  [15] 

0.94  [20] 

The  energy  balance  in  the  fuel  channel,  wall  1,  and  wall 
2  is  expressed  as 


3 

dx 


yj  M /-fcha-in  C i  Tf-fcha  O) 


-^SO^F-SRe{  Tf-FCHA  (*)  -  7wall2W}  =  0, 


(13) 


^F-So{7f-SoM  -  ^WalllWl  +  ^f-sre{Tsre(*) 


—  ^Walll  (■*)}  +  ^Wall^ALM 

and 


92  Twain  U) 

dx2 


=  0, 


(14) 


^F-SRe{TsRE  (x)  -  Twall2  (*)}  +  ^F-SRe{  Tf-FCHA  (*) 


~  Twall2  W)  +  ^Wall^-ALM 


a2  Twall2  (■%) 
dx2 


0, 


(15) 


where  the  heat  transfer  coefficient  in  the  fuel  channel  and 
the  stack  reformer,  /zf-sre,  and  the  heat  transfer  coefficient 
at  the  anode  side  of  the  tubular  SOFC,  /ip-so,  are  [16] 

ApA^uF 


^f-sre  = 


2/fcha  ’ 


and 


/if-so  = 


X^NuuP 

Vso 


(16) 


(17) 


The  molar  flow  rate  in  the  stack  reformer,  Mch4-sre,  varies 
according  to  [15] 


aMcH4-SRE(v)  (  -SSR  \ 

—  c - I  -  I 


3x 


=  ^SRE^SR  exp 


\RTske(x)J 


X  pch4-sre(v) 


0.9845 


(18) 


The  equilibrium  constant  of  the  shift  reaction,  ATshif,  is 
expressed  by 


TfsfflF{TsRE(*)}  = 


PCQ2  -  SRE  (*)  PR2  -SRE  (*) 
Pco-sre(-t)/?h2o-sre(v)  ’ 


(19) 


using  the  partial  pressure  of  each  gas  component  in  the  stack 
reformer.  The  mass  and  energy  balance  in  the  stack  reformer 
is  expressed  by 


3Mh2o-sre(v)  3MCh4-sreW  3MCq2-sre(v)  _  Q 


3a 


3a 


3a 


(20) 


9Mh2-SRe(v)  2  ^4Tch4-SRe(v)  _  3Mco2-SRe(v)  _  Q 
and  3a  3a  3a 

(21) 

-(Tf-pre-in  -  Tf-pre-out)^  Q Mi-pre-in 

1  a  _  ^CO-SReM  (  3MCh4-sreW  (  3MCo2-sre(v)  _  n 

+  (4TCh4-pre-in  -  ch4 -pre-out)  A Hsr  3a  +  3a  +  3a  _  U’ 

-  (Mco2-PRE-IN  -  4fc02-PRE-0UT)A/fsHIF  =  0.  (12)  (22) 
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and 


|r.  ^-sre  (x)  C  j  7s  re  (a) 


d2TSRE(x) 
dx 2 

9Mz-_sreW  1  , 


-  {TsreM  -  Th} 
9^ch4-sreW 


dx 


AHSr 


dMp  Q2-SRe(*) 
3a 


A//shif  -  <pso^f-sre{Tsre(v) 


-  Twain  (*)}  -  ^SO^F-SRe{TsreW  -  T\Vall2 CO }  =  0. 


(23) 


The  equilibrium  constant  of  the  oxidation  of  hydrogen, 
ATh2-oxe  that  of  carbon  monoxide,  A+o-oxe  and  that  of 
methane,  ATch4-oxi,  are  expressed  by 


TTh2-oxi{Tf-so(+)}  = 

Tfco-oxi{TF-soW}  = 
and 

TfcH4-OXI  {  Tf-so  C*0  }  = 


Ph2o-soW 


PU2-SO(x){po2-ANO(x)}0-5  ’ 

Pco2-so(x) _ 

PCO-SO (x){po2-ANo(x)}0  5  ’ 

{^H2o-soW}2i9co2-so(^) 

7>CH4-SO  (x)  {po2-ANO  (x)  }2 


(24) 

(25) 

(26) 


using  the  partial  pressure  of  each  gas  component  at  the  anode 
side  of  the  tubular  SOFC.  The  mass  and  energy  balance  at 
the  anode  side  of  the  tubular  SOFC  is  expressed  by 

dMn2o-so(x)  2 3MCH4-SoC* )  dMH2-So(*)  _  q 

dx  dx  dx  ~  ’  {  ) 

dMco2-so(x)  3Mch4-soW  ,  3Mco-soW 


dx 


+ 


dx 


+ 


dx 


=  0,  (28) 


and 


|  ^  3Tz-so  WQTp-so  W  J 


+  fCELL-SoWV  Q - 1  S°^  ^  -  ^S0^F-S0{Tf-S0W 

^  dx 

i 

-  Twain  Ml  -  <pso^f-so{Tf-soM  -  Tcell-soC*)}  =  0. 

(29) 

The  mass  and  energy  balance  at  the  cathode  side  of  the 
tubular  SOFC  is  expressed  by 

3MCh4-soM  .  3Mh2-soW  .  ^  cdMCo-so(x) 


+  0.5- 
3a 

dMp2- CATHW 
dx 


dx 


+  0.5- 


dx 


(30) 


and 

d 

dx 


{ 4T 02  -CATH  (x)  Cp2  Ta-SO  (*) 


+  4Tn2-AFM-InC'n2Ta-SoW}  +  Q)2Tceff-so(+) 


dMp2- CATHW 
dx 


^AFM^A-AFM-OUTS  {  T+-SO  W 


-  TafmM)  -  ^so^a-cath{ Ta-so  W 

-  Tceff-soM}  =  0, 


(31) 


where  the  heat  transfer  coefficient  between  the  air  supplied 
to  the  cathode  of  the  tubular  SOFC  and  the  cathode  of  the 
tubular  SOFC,  /ia-cath,  is  given  by 
^aA^Ua-cath<5so 


^A-CATH  = 


Y2  _  ,n2 

°SO  PALM 


(32) 


and  the  heat  transfer  coefficient  between  the  air  supplied 
to  the  cathode  of  the  tubular  SOFC  and  the  alumina  tube, 
/ia-aem-outs,  is  given  by 

^  A  Mi  A- AFM-OUTS  ^ALM 


^A-ALM-OUTS  = 


5SO 


^ALM 


(33) 


The  energy  balance  in  the  tubular  SOFC  is  given  by 

<PS0^F-S0{Tf-S0(v)  -  TcELL-SoWI  +  7r<$S0^A-CATH 
x  {  Ta-SO  W  -  TcELL-SoWI  -  {TcELL-SoW  -  Th} 

E  dMi- so  CO 

Ci - — - {TcELL-SoW  -  Th} 


X  Co 


dx 

3Mo2-cath(v)  3Mch4-soW 


3a 


+ 


3a 


■  A//ch4-oxi 


.  3Mh2-soW  a  .  3MCo-so(v) 

H - - - A  Hu2  -oxi  H - - - A  TTco-oxi 

3a  3a 

+  ^soOTano^ano  +  <Tcath^cath  +  <Tele^ele) 


32Tcell-so(v) 

3a2 


x(psoJso(x)Vsp. 


(34) 


The  current  density  of  tubular  SOFC,  Jso  CO*  is  given  [17] 
AF  3Mo2-cath(v) 


/so  CO  =  -- 


(35) 


7T(psP  dx 

The  voltage  of  the  tubular  SOFC,  +so>  is  given  by 
Vso  =  ESp(x)  -  0so  CO  -  Jso(x)Rc,  (36) 

where  the  electromotive  force  of  tubular  SOFC,  CsoCO*  is 
calculated  by  the  following  Nernst  equation  [9] : 


Eso(x)  =  - 


^Tcell-soC*) 
~AF 


In 


Po2-ano(x) 

/?o2-cathCO 


(37) 


The  ohmic  overpotential  of  tubular  SOFC,  0 so  00,  is  given 


by 

0soCO  =  ^{Tcell-soWI^soW, 


(38) 


where  the  total  resistance  of  the  tubular  SOFC, 
Rq { Tcell-S0 CO } >  which  is  the  summation  of  the  anode, 
cathode,  and  electrolyte  resistances,  is  expressed  by 

Ci 


^{Tcell-soOO}  =  7?Bexp 


Tcell-soCO 


+  C2 


(39) 
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where  the  values  for  C\  and  C2  were  estimated  to  be  5822 
and  —14.4,  respectively.  These  values  were  estimated  from 
the  I-V  characteristics  of  the  tubular  SOFC  made  by  SWPC 
at  1173  and  1273  K  [18],  which  had  no  voltage  drop  caused 
by  the  contact  resistance.  We  assumed  that  voltage  drop  in 
the  tubular  SOFC  made  by  SWPC  consists  of  only  ohmic 
overpotential  when  we  estimated  the  values  for  C\  and  C2. 

The  energy  balance  of  the  air  flow  inside  the  alumina  tube 
and  that  of  the  alumina  tube  are  expressed  by 
9 

-  —  {^A-SO-InC'mrTa-ALmW}  -  ^ALM^A-ALM-INS 
OX 

x  {Ta-almM  -  TalmOO}  =  0,  (40) 


^ALM^ A-ALM-OUTS  { Cx-SO  CO  -  TalmC*)} 

+  ^alm^a-alm-ins{^a-almW  -  7alm(*)} 


+  ^ALM^ALM^ALM 


d  TALM(x) 
dx2 


=  0. 


(41) 


The  heat  transfer  coefficient  inside  the  alumina  tube, 
/ia-alm-ins,  is  given  by 

^aNUa-alm-ins 


^A-ALM-INS  = 


<$ALM 


(42) 


The  energy  balance  in  the  combustor  is  expressed  by 
-A/7ch4-OXi(4^CH4-SO-OUT  +  4fcH4-PEstack-OUl7  AO 

-  A//h2-OXi(4/h2-SO-OUT  +  4fH2-PEstack-OUT/A0 

—  A/7coA^co-so-out  —  Q/N 

=  CaMa-ALM-In(7a-INS-ALM-COM-OUT  —  Tx-ALM-In) 

+  Co2Mo2-SO-OUt(7cEG-COM  -  ^a-so-out) 

+  Cn2Mn2-ALM-In(^CEG-COM  —  ^a-so-out) 

+  ^  QM.-SO-out(7ceg-com  -  7f-so-out) 

i 

+F,  Ci  Afj-PEstack-OUT  (  ^CEG-COM 
i 

-  ^F-PEstack— OUT ) / N,  (43) 

where  the  air  temperature  inside  the  alumina  tube  at  outlet  of 
the  combustor,  7a-ins-alm-com-out,  is  calculated  by  [19] 

C'a4^a-alm-in(^a-ins-alm-com-out  —  ^a-alm-in) 

=  ^COM^COM 

(Tceg-com  -  7a-ALM-In) 

—  (7a-so-out  —  Tx-ins-alm-com-out) 


ln{(TcEG-COM  -  7a-ALM-In)/ 

(^A-SO-OUT  -  7a-INS-ALM-COM-OUt)} 


(44) 


The  molar  flow  rate  of  the  combustion  exhaust  gas,  Mqeg, 
is  given  by 


The  dc  output  of  the  tubular  SOFC,  Wso,  is  given  by 

/»/cell 

Wso  =  ncpsoVso  /  Jso(x)dx.  (46) 

Jo 

The  molar  flow  rate  of  the  anode  exhaust  gas  for  the  PEFC 
stack,  Mape,  is  given  by 

Mape  =  ^aeg-so  -  Ma.  (47) 

The  molar  flow  rate  of  the  methane  fed  to  the  SOFC  stack, 
A^CH4-S0stack-IN>  the  molar  flow  rate  of  the  air  supplied  to 
the  SOFC  stack,  MA-S0stack-IN>  the  dc  output  of  the  SOFC 
stack,  Wsostack-dc » the  molar  flow  rate  of  the  anode  exhaust 
gas  for  PEFC  stack  from  the  SOFC  stack,  MAPE-S0stack, 
and  the  molar  flow  rate  of  the  combustion  exhaust  gas  from 
the  SOFC  stack,  McEG-S0stack>  are  given  by  the  following 


equations: 

A^CH4-SOstack-IN  =  4fcH4 -PRE-IN  N,  (48) 

A^X-SOstack-IN  =  4fA-ALM-INW  (49) 

WsOstack-dc  =  WsoN,  (50) 

A^APE-SOstack  =  MaPeN,  (51) 

and 

A^CEG-SOstack  =  MqegN.  (52) 


3.2.  SOFC-PEFC  system 

The  dc  output  of  the  PEFC  stack,  WpEstack-dc>  is  calculated 
from  the  fuel  utilization  rate  of  the  PEFC  stack,  t/fPEstack> 
the  voltage  of  the  PEFC,  Vpe,  and  the  molar  flow  rate  of  the 
hydrogen  in  the  exhaust  gas  from  the  CO  selective  oxidizer, 
A^h2 -PEstack-iN >  using  the  following  equation: 

TI7  2Mn2-PEstack-IN  f/fPEstack  VpE  F 

WpEstack-dc  =  - — - ,  (53) 

where  MH2-PEstack-iN  is  given  by 


4^H2  -PEstack-IN  =  4fR2-APE-SOstack  +  0.99McO-APE-SOstack- 

(54) 

The  temperature  of  the  air  supplied  to  the  SOFC  stack, 
^A-S0stack-iN » is  calculated  based  on  the  following  equations: 

Ca A-SOstack-IN ( ^A-SOstack-IN  “  ?AMb) 

=  CcOM^CEG-SOstack(7cEG-SO-OUT  “  TcEG-HE-OUt)- 

(55) 
and 


MqeG  =  (4/h20-S0-0UT  +  4fH2O-PEstack-OUT/A0 
+  (Mco2-SO-OUT  +  4fcO2-PEstack-OUT/A0 
+  (Mo2-SO-OUT  +  ATn2-S0-0UT  ) 
-0.5(Mh2-SO-OTU  +  4fH2-PEstack-OUT/A0 
— 0.5Mco-so-out- 


Ca  A-SOstack-IN  ( ^A-SOstack-IN  —  ?AMb) 

(^CEG-SO-OUT  -  ^A-SOstack-IN ) 

—  (7ceg-he-out  —  Txmb) 

=  a  HE  4  HE - 

ln{(7cEG-SO-OUT  -  ?A-SOstack-IN )  / 
(7ceg-he-out  -  ?amb)} 


(56) 


(45) 
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The  ac  output  of  the  SOFC  stack,  Wsostack-ac,  and  that  of 
the  PEFC  stack,  WpEstack-ac,  are  given  by 

fksOstack-ac  =  V^(l  —  ^  oAux  ~  SOstack)lFso  stack- dc>  (57) 

and 

WpEstack-ac  =  ^(1  -  T0AUX-PEstack)  WpEstack-dc,  (58) 

where  the  ratio  of  the  power  consumption  in  the  auxil¬ 
iary  machine  to  the  gross  ac  output  of  the  SOFC  stack, 
T'oAUX-sostack » is  calculated  by 

T  __  _  ^^oblow  I  r 

^c>Aux-SOstack  —  TT  i 


u, 


OQTHER • 


(59) 


°SOstack 


Oxygen  utilization  rate  in  SOFC  stack  [%] 


The  ac  output  of  the  SOFC-PEFC  system,  Wsys  ,  and  the 
electrical  efficiency,  77,  are  calculated  by 

Wsys  =  ^sOstack-ac  +  WpEstack-ac »  (60) 

and 

100(WpEstack-ac  +  WsOstack-ac)  raw 

rj  = - .  (61) 

^CH4-SOstack-IN  A  #ch4-OXI 

4.  Results  and  discussion 

4.1.  Effectiveness  of  simulation 

The  ac  output  and  the  electrical  efficiency  of  the  SOFC 
cogeneration  system  made  by  SWPC  were  simulated  by  us¬ 
ing  the  simulation  model  shown  in  Fig.  1.  We  compared  the 
simulation  result  with  the  actual  performance  demonstrated 
by  SWPC  [2]  (see  Table  2).  The  former  is  in  good  agree¬ 
ment  with  the  later.  We  can  conclude  that  our  simulation  is 
effective  at  estimating  ac  output  and  electrical  efficiency. 

4.2.  Effect  of  oxygen  utilization  rate  in  the  SOFC  stack 

The  dependence  of  the  ac  outputs  of  the  SOFC  stack, 
the  PEFC  stack,  and  the  SOFC-PEFC  system  on  the  oxy¬ 
gen  utilization  rate  in  the  SOFC  stack  is  shown  in  Fig.  3. 
The  maximum  temperature  in  the  SOFC  stack  is  less  than 
1273  K  when  the  oxygen  utilization  rate  in  the  SOFC  stack 
is  less  than  30%,  because  a  large  amount  of  exhaust  heat  of 
the  SOFC  stack  is  discharged  from  the  stack  when  a  large 
amount  of  air  supplied  to  the  SOFC  stack.  The  methane 
with  a  molar  flow  rate  exceeding  the  molar  flow  rate  needed 
to  keep  the  average  current  density  of  the  tubular  SOFC 
at  2000  A/m2  is  reformed  in  the  stack  reformer  while  the 
maximum  temperature  in  the  SOFC  stack  is  maintained  at 


Table  2 

Comparison  of  ac  output  and  electrical  efficiency 


ac  output  (kW) 

Electrical  efficiency(%) 

Simulation 

114 

47 

SWPC  [2] 

109 

46 

Fig.  3.  Dependence  of  ac  outputs  of  the  SOFC  stack,  PEFC  stack,  and 
SOFC-PEFC  system  on  the  oxygen  utilization  rate  in  the  SOFC  stack.  Tm : 
maximum  temperature  in  SOFC  stack,  J so- awe'-  average  current  density  in 
SOFC  stack,  Vfe:  cell  voltage  in  PEFC  stack  and  £/fPEstack :  fuel  utilization 
rate  in  PEFC  stack. 


1273  K  when  the  oxygen  utilization  in  the  SOFC  stack  is 
over  30%.  This  is  because  the  amount  of  exhaust  heat  of  the 
SOFC  stack  discharged  with  the  air  decreases  with  increas¬ 
ing  oxygen  utilization  rate  in  the  SOFC  stack  and  because 
that  used  effectively  as  the  reaction  heat  of  the  steam  reform¬ 
ing  reaction  of  methane  in  the  stack  reformer  increases.  The 
ac  output  of  the  PEFC  stack  increases  with  increasing  oxy¬ 
gen  utilization  rate  in  the  SOFC  stack  because  the  methane 
used  for  the  power  generation  of  the  PEFC  stack  increases. 
The  ac  output  of  the  SOFC  stack  slightly  increases  with  in¬ 
creasing  oxygen  utilization  rate  in  the  SOFC  stack.  This  is 
because  the  cell  voltage  of  the  tubular  SOFC  increases  since 
the  partial  pressures  of  hydrogen  and  carbon  monoxide  in 
the  steam-reformed  gas  at  the  anode  of  the  tubular  SOFC  in¬ 
crease,  and  because  the  power  consumption  of  the  air  blower 
decreases.  Consequently,  the  ac  output  of  the  SOFC-PEFC 
system,  which  is  the  sum  of  the  ac  outputs  of  SOFC  and 
PEFC  stacks,  increases  with  increasing  oxygen  utilization 
rate  in  the  SOFC  stack.  The  amount  of  steam  in  the  recycled 
gas  for  the  steam  reforming  is  insufficient  and  the  S/C  ratio 
in  the  mixed  gas  fed  to  the  pre-reformer  is  less  than  3.0  when 
the  oxygen  utilization  rate  in  the  SOFC  stack  is  more  than 
70%.  The  ac  output  of  the  SOFC-PEFC  system  is  175  kW 
when  the  oxygen  utilization  rate  in  the  SOFC  stack  is  70%. 

The  dependence  of  the  electrical  efficiency  of  the 
SOFC-PEFC  system  on  the  oxygen  utilization  rate  in  the 
SOFC  stack  is  shown  in  Fig.  4.  The  electrical  efficiency 
of  the  SOFC-PEFC  system  increases  with  increasing  oxy¬ 
gen  utilization  rate  in  the  SOFC  stack.  This  is  because  the 
amount  of  exhaust  heat  of  the  SOFC  stack  used  as  the  reac¬ 
tion  heat  of  the  steam  reforming  reaction  increases  and  be¬ 
cause  that  discharged  from  the  stack  with  the  air  decreases. 
The  electrical  efficiency  of  the  SOFC-PEFC  system  is  54% 
when  the  oxygen  utilization  rate  in  the  SOFC  stack  is  70%. 
This  value  is  far  higher  than  the  electrical  efficiency  at 
109  kW  ac  of  the  SOFC  cogeneration  system  demonstrated 
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Oxygen  utilization  rate  in  SOFC  stack  [%] 


Fig.  4.  Dependence  of  the  electrical  efficiency  of  the  SOFC-PEFC  system 
on  the  oxygen  utilization  rate  in  SOFC  stack.  Tm :  maximum  temperature 
in  SOFC  stack,  /so- aye-  average  current  density  in  SOFC  stack,  Vpe:  cell 
voltage  in  PEFC  stack  and  t/fPEstack :  fuel  utilization  rate  in  PEFC  stack. 


Fig.  6.  Dependence  of  the  electrical  efficiency  and  ac  output  of  the 
SOFC-PEFC  system  on  cell  voltage  in  the  PEFC  stack.  Tm:  maximum 
temperature  in  SOFC  stack,  /so-aye-  average  current  density  in  SOFC 
stack  and  f/fPEstack-  fuel  utilization  rate  in  PEFC  stack. 


by  SWPC,  46%  [2],  but  lower  than  that  at  220  kW  ac  of  the 
SOFC-GT  system  designed  by  SWPC,  57%  [4]. 

4.3.  Effect  of  fuel  utilization  rate  and  cell  voltage  in  PEFC 
stack 

The  fuel  utilization  rate  of  the  PEFC  stack  in  a  PEFC  sys¬ 
tem  is  low  because  the  large  amount  of  hydrogen  in  anode 
exhaust  gas  has  to  be  burnt  to  generate  the  heat  needed  for 
the  steam  reforming  reaction  of  methane  and  steam  genera¬ 
tion  [22].  On  the  other  hand,  the  fuel  utilization  rate  of  the 
PEFC  stack  in  the  SOFC-PEFC  system  can  be  increased 
compared  with  the  PEFC  system,  because  the  temperature 
of  the  exhaust  heat  from  the  SOFC  stack  is  about  1273  K 
and  this  heat  can  be  used  for  the  steam  reforming  reaction  of 
methane  and  steam  generation.  The  dependence  of  the  elec¬ 
trical  efficiency  and  the  ac  output  of  the  SOFC-PEFC  sys¬ 
tem  on  the  fuel  utilization  rate  in  the  PEFC  stack  is  shown 
in  Fig.  5.  The  electrical  efficiency  and  the  ac  output  of  the 


SOFC-PEFC  system  increase  with  increasing  fuel  utiliza¬ 
tion  rate  in  the  PEFC  stack.  The  electrical  efficiency  and  the 
ac  output  of  the  SOFC-PEFC  system  are  56%  and  180kW, 
respectively,  when  the  fuel  utilization  rate  in  the  PEFC  stack 
is  85%.  The  electrical  efficiency  of  56%  is  almost  equal 
to  that  at  220  kW  ac  of  the  SOFC-GT  system  designed  by 
SWPC,  which  is  57%. 

The  dependence  of  the  electrical  efficiency  and  the  ac 
output  of  the  SOFC-PEFC  system  on  the  cell  voltage  in 
the  PEFC  stack  is  shown  in  Fig.  6.  The  electrical  efficiency 
and  ac  output  of  the  SOFC-PEFC  system  increase  with  the 
cell  voltage  in  the  PEFC  stack.  The  electrical  efficiency  and 
ac  output  of  the  SOFC-PEFC  system  are  59%  and  190  kW, 
respectively,  when  the  cell  voltage  in  the  PEFC  stack  is 
0.85  V.  The  electrical  efficiency  of  59%  is  equal  to  that  at 
1014  kW  ac  of  the  SOFC-GT  system  designed  by  SWPC 
[4]. 
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Fig.  5.  Dependence  of  the  electrical  efficiency  and  ac  output  of  the 
SOF-PEFC  system  on  the  fuel  utilization  rate  in  PEFC  stack.  Tm:  max¬ 
imum  temperature  in  SOFC  stack,  /so-ave-  average  current  density  in 
SOFC  stack  and  Vpe:  cell  voltage  in  PEFC  stack. 


5.  Summary 

We  evaluated  the  200  kW-class  SOFC-PEFC  system  by 
numerical  simulation.  The  main  results  are  as  follows. 

1.  The  electrical  efficiency  of  the  SOFC-PEFC  system  in¬ 
creases  with  increasing  oxygen  utilization  rate  in  the 
SOFC  stack.  This  is  because  the  amount  of  exhaust  heat 
of  the  SOFC  stack  discharged  from  the  SOFC  stack  with 
the  supplied  air  decreases  as  the  oxygen  utilization  rate 
in  the  SOFC  stack  increases  and  because  that  used  effec¬ 
tively  as  the  reaction  heat  of  the  steam  reforming  reaction 
of  methane  in  the  stack  reformer  increases. 

2.  The  electrical  efficiency  at  190  kW  ac  of  the  SOFC-PEFC 
system  is  59%  when  the  fuel  utilization  rate  and  the  cell 
voltage  in  the  PEFC  stack  are  85%  and  0.85  V,  respec¬ 
tively.  This  electrical  efficiency  is  equal  to  that  of  the 
SOFC-GT  system  at  1014  kW  ac. 
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